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ABSTRACT 

We report the results of a multi-waveband analysis of the masses and luminosities of 
~600 galaxy groups and clusters identified in the maxBCG catalogue. These data are 
intended to form the basis of future work on the formation of the “mi 2 gap” in galaxy 
groups and clusters. We use SDSS spectroscopy and g, r and i band photometry to 
estimate galaxy group/cluster virial radii, masses and total luminosities. In order to 
establish the robustness of our results, we compare them with literature studies that 
utilize a variety of mass determinations techniques (dynamical, X-ray, weak lensing) 
and total luminosities estimated in the B , r, i, and K wavebands. We also compare 
our results to predictions derived from the Millennium Simulation. We find that, once 
selection effects are properly accounted for, excellent agreement exists between our 
results and the literature with the exception of a single observational study. We also 
find that the Millennium Simulation does an excellent job of predicting the effects 
of our selection criteria. Our results show that, over the mass range ~ 10 13 — 10 15 
Mq , variations in the slope of the mass-luminosity scaling relation with mass detected 
in this and many other literature studies is in part the result of selection effects. 
We show that this can have serious ramifications on attempts to determine how the 
mass-to-light ratio of galaxy groups and cluster varies with mass. 
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1 INTRODUCTION 


In 1933 Fritz Zwicky applied the virial theo rem to a galax y 
cluster (the Coma cluster) for the first time (j Zwic .kv 1 19331) . 
He concluded that, in order to explain the cluster’s dynam¬ 
ics, the average density of the cluster had to be several hun¬ 
dred times greater than that indicated by estimates of the 
mass of luminous material observed in the cl uster’s galaxie s 
alone. Confirmed later in the Virgo cluster (Smith Il93a) . 
these were, of course, the first direct observations indicating 
the presence of dark matter (Zwicky’s “Dunkle Materie”). 


Zwicky’s methodology of measuring the mass of clus¬ 
ters and comparing this mass to the total light emitted by 
the galaxies in the cluster is still in use today. Indeed, the 


measurement of mass-to-light ratios in galaxy groups and 
clusters has become commonplace since these discoveries, as 
knowledge of how the mass of baryonic and non-baryonic 
matter are distributed within clusters provide important 
clues as to how they, and the galaxies within them, were 
formed. 

However, the techniques and data involved in measur¬ 
ing both the masses and luminosities of clusters have become 
significantly more sophisticated since Zwicky’s original dis¬ 
covery. For instance, with the coming of large telescopes 
and large scale surveys, the measurement of group/cluster 
luminosities has advanced, permitting direct observation of 
galaxies further down the luminosity function and to higher 
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redshifts, while mass estimates can now be based upon 
cluster dynamics (as performed by Zwicky and Smith), or 
in more recent advances, X-ray properties of cluster halos 
or gravitational lensing. However, due to the variations in 
methodologies adopted in works in the literature, the re¬ 
ceived wisdom has become that it is difficult to make direct 
comparisons between studies. 

In this work we examine this issue, comparing the re¬ 
sults of our analysis of the masses and luminosities of ~600 
galaxy groups and clusters identified in the maxBCG cata¬ 
logue with the results of the Millennium Simulation and five 
other recent studies in the literature that use a variety of 
methodologies and data sources. Our analysis is performed 
in four wavebands ( B , r, i and K) and includes systems 
of mass ranging from ^10 13 Mq to 10 15 Mq. We examine 
both the extensive areas of agreement, as well as areas of 
disagreement, between these studies and identify the issues 
that must be considered before making comparisons between 
studies. 

In future papers we will use the data presented here to 
perform an analysis aimed at identifying the group/cluster 
properties driving differences in the “mi 2 gap” (the lumi¬ 
nosity difference between the brightest and second brightest 
galaxies within the central regions of groups and clusters). 

Throughout this paper both our data and literature 
data are presented assuming a Ho=70 km s -1 Mpc -1 , 
Da=0.7, Dm=0.3 cosmology. For expressing luminosities in 
solar units we used the values M r ,0=4.67 mag, M; i q=4.48 
mag, Mb,0=5.33 mag, Mif : 0=3.28 mag. 


clusters studied have minimal impact on distance determi¬ 
nations and their associated distance moduli, and thus en¬ 
suring accurate luminosity estimates. A minimum of 6 spec¬ 
troscopically confirmed members was required as the mass 
estimates made in this work are based on velocity dispersion 
estimates whose errors are large when the number of spec¬ 
troscopically confirmed members is low. Indeed, even with 
a sample of 6 spectroscopically confirmed members, errors 
in log(mass) estimates are of order of 0.4 dex (or a factor of 
2.5 in mass). 


2.1 Data from the literature 


In this work we compare our results to those o f five recent 
studies . These studies ar e: |GjrmTli_gt_ahJ fl2002|), Ekeet aTTI 
(l2004bh. Ramell a et al. I (12004 1. Popesso et al. I (2007) and 
ISheldon et al. I ( 2009fl . These studies span four wavebands 
( B , r, i and K). In each case the results of these studies 
were, if necessary, converted to the Ho=70 km s -1 Mpc -1 , 
Da=0.7, Dm =0.3 cosmology used throughout this paper. 
Where necessary, small corrections were also made to ensure 
that the same solar luminosities were used in the estima¬ 
tion of group/cluster luminosities and that the data are k- 
corrected to the same redshift (z = 0). We note that the lit¬ 
erature studies above are generally based on groups/clusters 
at lower redshift than our data, and a minimum of 10 spec¬ 
troscopically confirmed members is generally applied (rather 
than the 6 used here). The studies using dynamical mass es¬ 
timates therefore generally possess average mass errors that 
are smaller than ours by, in the worst case, ~0.1 dex. 


2 DATA 

Our sample was selected from grou ps identified in the 
maxBCG catalogue of iKoester et al.~l (|2007bh . This cata¬ 
logue was constructed from the SDSS photometric survey 
using a cluster finding algorithm based on three well de¬ 
fined properties of galaxies in clusters: spatial clustering, 
the presence of a “red-sequence” and the presence of a 
BCG (brightest cluste r galaxy) at the cente r of the cluster 
I Koester et al. If2007al l. Using this method IKoester et al. I 
I 2007bl) identi fy more than 13,0 00 clusters. Using DR9 of 
the SDSS-III ( Ahn et ah~l 120121 1 , we then selected all the 
clusters from the maxBCG catalogue that have a BCG with 
a spectroscopic redshift in the redshift interval z = 0.05 to 
2 = 0.16 and possessing 6 or more spectroscopically con¬ 
firmed members (as described in Section [Oil . 

The upper end of the selected redshift range was cho¬ 
sen in order to ensure a substantive sample at all values of 
the “mi 2 gap”, as it is necessary to probe out to redshifts 
of ~0.15 in order to obtain a reasonably large sample of 
the rare “fossil groups” (defined as systems with 77112 > 2). 
The lower end of the redshift range was chosen to minimize 
the impact of the redshift dependent selection criteria of the 
SDSS spectroscopic survey which only sampled galaxies with 
apparent magnitude brighter than 17.7 mag in the r band, 
and to ensure that the peculiar velocities of the groups and 


2.2 The Millennium Simulation 

To supplement our analysis and provide a means of test¬ 
ing the impact on observations of errors and selection ef¬ 
fects, we also use the galaxy recession velocities (for dy¬ 
namical mass estimation) and luminosity data from the 
Millennium Run dark matter simulation in the WMAP1 
cosmology (Ho=73 km s -1 Mpc -1 , Da=0.75, Dm=0.25) 
from Springel et al. (2005), together with the semi-analytic 
galaxy modeling technique as described in Guo et al. 
(2011). The semi-analytic galaxy catalogs were converted 
into an obs erver’s light-cone geomet ry matched to the SDSS 
survey by iHenriques e t al. I (l2012fl using the model s of 
iBruzual fe Chariot I ll2003lf. assuming a Chabrier ]Hqq 3) ini- 
tial mass function. IHenriques et al. I (120121 1 estimated total 
luminosities in a broad range of photometric bands from 
the B band to the infrared. In this work we use the r band 
data only. These data were converted to the Ho=70 km s -1 
Mpc -1 cosmology used throughout this paper. The correc¬ 
tions for Da and Dm were deemed insignificant for this work 
as they cause differences of only ~0.01 mag and ~0.4% in the 
distance moduli and angular sizes, respectively. Data for all 
galaxies brighter than 22 mag in the r band (the photomet¬ 
ric limit of the SDSS data used in our analysis) in simulated 
halos at redshifts less than 0.16 were obtained, and the po- 
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sitions and r band luminosities of each galaxy, along with 
the M 200 of the parent halos, were tabulated. 


2.3 Spectroscopic selection and analysis 

The identification of spectroscopically confirmed cluster 
members was performed by an iterative procedure. First, 
each cluster was assumed to possess a velocity dispersion 
(a) of 600 km s -1 . R200 was the n estimated assuming the 
relation of lCarlberg et al.l ljl997h : 


R200 = 


\/3 a 
10 H{z) 


Mpc. 


( 1 ) 


We then searched the SDSS spectroscopic catalogue for 
all galaxies within the estimated R200 that possess redshifts 
differing from the BCG redshift by less than 1500 km s _1 
(i.e. ±2.5 <t). The value 2.5cr was selected in order to allow 
for the inherent scatter expected from the (almost Gaussian) 
distribution of velocities observed in clusters, while at the 
same time minimizing the likelihood of selecting interlopers. 
Once all the members within R200 were identified their red- 
shi fts were used to calculate a new velocity dispersion using 
the lBeers et al. IJ 1 99(1 ) bi-weigh t estimator fo r syste ms with 
more than 10 members, and the lBeers et al. I lll990li gapper 
method for systems with 10 or less members. For cases where 
the initial velocity dispersions were an underestimate this re¬ 
sulted in an increase in the velocity dispersion (and hence 
R200) estimate, as well as increasing the permitted velocity 
range (±2.5<r n ew). Conversely, for systems in which the ini¬ 
tial values were an overestimate the new values were lower. 
This process was iterated a number of times. We found that, 
generally, convergence was achieved after only two or three 
iterations. However a total of 10 iterations were performed. 
This permitted the use of the value from the last (10th) it¬ 
eration as the final result for each group/cluster, while the 
RMS scatter in the last 6 iterations was taken as an addi¬ 
tional uncertainty in the velocity dispersion and which was 
added in quadrature to the formal velocity dispersion error. 
We note that this additional error is zero for the majority 
of systems in which the procedure converged. 

At this stage, groups/clusters with less than 6 spectro¬ 
scopically confirmed members within R200 were discarded 
from the sample. Clusters for which the average recession 
velocity of all the galaxies except the BCG differed from the 
recession velocity of the BCG by more than 400 km s -1 
were also excluded, as these have a high probability of be¬ 
ing either highly disturbed (i.e. non-virialized) systems or 
false/contaminated detections due to line-of-sight effects. 


2.4 Mass estimation 

The masses of the remaining maxBCG clusters were esti¬ 
mated using: 

3 2 

A /200 = ^ 0200 ^ 200 - (2) 


This equation can be reformulated using Equation [Q to 
give a more convenient form: 

log(M 2 oo) = 31og(cr2oo) + 6.24 — log(E(z)), (3) 

where, for a flat Universe: 

E{z) = V^m(1 + z) 3 + Qa. (4) 

However, it is helpful to note that log(E(z)) is small 
(~0.02 at the z = 0.1 median of the data presented here) 
and varies by only ±0.02 over the whole redshift range of 
the clusters reported here. The resultant masses spanned 
the range from ~10 13 to ~10 15 Mg. 


We note that four of the five literature studies to which 


we compare our results 1 

Girardi et al. 200^ Eke et al. 

2004b: Ramella et al. 2004 

: Popesso et al. 200% use Equa- 


tion \2\ to estimate masses. However, both Girardi et al. I 


d2002ll and lPopesso et al~l (l2007ll make corrections for a sur- 
face pressure t e rm. A comparison of the masses presented in 
iGirardi et al. I ((20021 4 with masses estimated using their ve¬ 
locity dispersions and Equation 0 yields an offset of 0.00 
dex with RMS scatter of only 0.12 dex. This demonstrates 
that this term i s small and of no significance to the results 
presented here. lEke et al~l (12004431 1 , on the other hand, use 
the spatial distribution of the galaxies in each system to 
estimate R 200 rather than Equation [T] basing their values 
of radius on a cali bration agai nst cosmological si mulations 
llEke et al. I l2()04al 'F1. Finally, dSheldon et al. 1 120091 1 utilize 
SDSS i band photometry and masses derived from weak 
gravitational lensing. 

As a first test of our dynamical mass estimates, we 
performed a search of the literature for groups/clusters 
in our study for which masses derived from X-ray prop¬ 
erties have either been reported or can be estimated. A 
total of 24 groups from our study were found to have 
masses based on a variety of X-ray properties reported 
in th e literature (|Dqnahue_eL_ai~| |2005|; Uia^_Chen_fc_ChenJ 


2006|_ Lemze_etah 200S Zibettj^jhemh^ PraTt_ 2009|; 


Akamatsu et al . I 1201 ll: IProctor et al. 201 ll : iMiller et al. I 


2012; ^wCTS_et_M. 2014). For the major ity of of these stud- 
Donahue et al. 20 1 Zibettij^ierini&Pratt^ 2009|; 


Aka matsu et al. I 120111 " IProctor et al. 201 it IMiller et al. I 


2012 ; Owers et al. Il20j~4h X-ray masses were estimated from 
X-ray temperatu re using the the Tier 1 relation from 
ISun et al. I (|2009|) (their Table 6 ) . How ever, the X-ray masses 
taken from I Jia, Chen fe Chen I (|2006i l were deived using X- 
ray temperat u re an d electron density, while in the case of 
iLemze et al. I (|2008l f . weak lensing and the X-ray luminos¬ 
ity profile were employed to determine cluster masses. The 
masses derived from all the above studies assume hydrostatic 
equilibrium. 


1 N.B. We have converted the bj given in I Eke et aT~l d2004bl 'l 
to the B band using the relation Lt, ;j = 1.1 Lg given by their 
Equation 4.3 
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The results of the comparison of dynamical masses to 
those derived from X-ray data are shown as solid points in 
Fig- EO Also shown in this figu re as op en sq uares are to¬ 


tal cluster masses derived from Zhang ct al 


sq uares c 
■ ( 20111 ) 


X-ray 


gas m asses assuming the relation derived by Mahdavi et al. I 
(|2013f l between X-ray gas masses and total cluster masses. 
This relation was calibrated by comparing gas masses from 
X-ray analyses to the results of total cluster masses from 
weak lensing. Dynamical masses for these d ata were de- 
rived using the velocity dispersions provided in lZhang et al. I 
(120111') and Equation [3] Velocity dispersions and masses 
derived from X-ray temp eratu re were also taken from 
lOsmond fe Ponman I d2004l) and IWu. Xue fe Fang~l dl999ll . 
For these data, dynamical masses were derived using Equa¬ 
tion 02 while X- ray masses w ere derived from the X-ray tem- 
peratures using SunjeLah. d2009l) as above. The re s ults o f 
lOsmond fe Ponman I ( 2004 1 and IWu. Xue fe Fang! dl999l l 
are shown in Fig. [T| as open circles and triangles respec¬ 
tively. Again, hydrostatic equilibrium is assumed in all cases, 
and all the above results are for masses within J? 200 - Where 
necessary, we assumed f?20o=l-54f?5oo, consis tent with the 
values of R200 and -R 500 for clusters given in lEttori et al. I 
d20ld') . yielding JV/ 2 oo = 1.46M 5 oo. 


Now, given the inhomogeneity of these data, and the 
need to convert between M500 and M200, it would be un¬ 
wise to over-interpret this plot. However, we note that in 
Fig.m systems with masses above ~10 14 Mq, although run¬ 
ning parallel to the one-to-one line, the X-ray masses ap¬ 
pear systematically offset from the dynamical masses by 
~0.15 dex. On the other hand, the systems with masses 
lower than ~10 14 Mq (which are from lOsmond fe Ponmanl 
1 2004i f 1. the data straddle the one-to-one line, but appear to 
exhibit a slope greater than one._This_change in gradient at 
low masses was noted in lOsmond fc Ponman I d2004ll when 
they considered their Lx — cr and a — Tx plots (their Sec¬ 
tions 7.2 and 7.3) and was also noted in lHelsdon fc Ponman I 
iteOOdf . A break in the slope of the X-ray scaling relations 
at a temperature of about 1 keV (corresponding to a mass 
of ^10 14 Mq) is also evident in a number o f studies (e.g. 
IXue fe Wu I (I200ol l and lHarrison et al~l ll2012lff. How ever, i t 
has been disputed in other works, e.g. iMulchaev I (2000). 
Therefore, there remains no real consensus as to whether 
the apparent change in gradient in the Lx — a and a — Tx 
scaling relations represents a true physical effect or is rather 
an artifact of large errors an d sel ection effects in low mass 
systems fsee lOsmond fe Ponmanl J2004l l for a more detailed 
description of the controversy). 

It is also interesting to perform a similar comparison us¬ 
ing the Millennium Simulation data, in this case comparing 
dynamical mass to the M 200 values given in the database. 
The masses presented in the database were calculated by 
simply adding the masses of all the particles within the re¬ 
gion of each halo where the average density is 200 times 
the critical density. In order to calculate dynamical masses 
we applied the same analysis as that applied to our ob¬ 
servational data. i.e. assuming that the brightest cluster 
galaxy lay at the center of the cluster and measuring re- 



Figure 1. Comparison of dynamical mass, estimated as detailed 
in the text, with masses derived from X-ray data. Black dots 
are clusters in our study (whose dynamical masses are estimated 
here) with X-ray masses from the literature. Open squares are de- 
r ived fr om the velo city dispersion and X-ray gas mass data from 
IZhang et al. I ll201lh (see text). Open circles and triangles are de- 
ri ved from the velocity dispe rsion and X-ray temperature data 
of I Osmond &; Ponman I i 20041) and lWu, Xue fc Fang I (119991 ). re¬ 
spectively. The dashed line shows the one-to-one locus. At masses 
greater than ~ 10 14 Mq there is a clear, constant offset between 
X-ray and dynamical masses, while below this mass there is a 
suggestion of a change in slope in the correlation. 


ces sion velocities with respect to the central galaxy using 
the Beers et al. I (|l990l ) bi-weight and gapper estimators as 
described in Section noito estimate velocity dispersion. Dy¬ 
namical masses are then estimated using Equation [3] This 
analysis was performed under the assumption of a variety 
of selection criteria, including those describing our observa¬ 
tional data, which we recall were selected to have redshifts 
between 0.05 and 0.16, and have 6 or more members with ap¬ 
parent magnitudes brighter than the 17.7 mag within i? 200 - 
For each analysis, the median of the dynamical mass data 
in bins of Millennium Simulation mass were calculated. We 
note that our use of the median in these analyses ensures 
that the effects of the asymmetrical errors in dynamical mass 
are essentially eliminated. 

The comparison between Millennium Simulation masses 
and dynamical masses is shown in Fig. [2] In this plot the 
dynamical masses are those for all systems with a minimum 
of 6 members, but with no redsliift or member luminos¬ 
ity selection limits applied. Reasonable agreement is found. 
At masses greater than ^10 14 Mq the Millennium Simu¬ 
lation masses show an offset from the median dynamical 
masses of ~0.1 dex, similar to the trend in the X-ray data 
in Fig. □ Below ~10 14 Mq, the correlation between Mms 
and Md S ti steepens, with the gradient increasing from 1.05 
above ~10 14 Mq to 1.17 below this mass. Although statis- 
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Figure 2. Comparison of masses using our dynamical mass esti¬ 
mation method (see text) with masses from the Millennium Sim¬ 
ulation. The thick red line shows the one-to-one locus, while the 
thin lines show the median and 15/85 percentiles of the dynami¬ 
cal masses in bins of 0.1 dex in Millennium Simulation mass. As 
well as an offset at high masses, the data exhibit and increase in 
slope at masses below 10 14 Mq. 

tically only marginally (~2cr) significant, this steepening at 
low masses is again at least qualitatively similar to the X-ray 
data. Some care must be taken in the interpretation of these 
differences, since given the large uncertainties in the semi- 
analytic models that describe the baryonic physics, there is 
no guarantee that the Millennium Simulation data match 
the real Universe in this regard. However, it is interesting 
to note that offset at high masses and the increased gra¬ 
dient at lower masses are at least qualitatively consistent 
with the trends in the X-ray data. This therefore raises the 
possibility that this represents a real systematic bias in dy¬ 
namical masses, perhaps due to a bias in the distribution of 
galaxies with respect to the underlying dark matter halo or 
perhaps representing a breakdown in the assumptions made 
in deriving dynamical masses (i.e. the assumption of full 
virialization and/or that the systems are well described as 
isothermal spheres). 

2.5 Photometric analysis 

Having established estimates of the masses of the clusters 
in our sample we must next estimate the total luminosity of 
each of the systems. Therefore, for each of the clusters se¬ 
lected above, the g, r and i band photometry of all galaxies 
within 4f?200 of the BCG were obtained from the SDSS pho¬ 
tometric survey. We use the SDSS model magnitudes which 
were estimated using a combination of De Vaucouleurs and 
exponential profiles. Galaxies with SDSS spectroscopic red- 
shifts were k-correct ed to zero redshift using th e colour de¬ 
pendent functions of IChilingarian et al. I J20ld ). 


For each cluster, and in each waveband ( g , r and i), 
the total luminosity of the galaxies within R200 was then 
estimated in a two step process. First, the luminosities of 
all the spectroscopically confirmed galaxies within -R 200 of 
each group were co-added to derive the total luminosity. 
Next, the total luminosity of all galaxies in the photomet¬ 
ric survey with apparent magnitudes fainter than 17.7 mag 
in the r band (the lower limit of the spectroscopic sur¬ 
vey) and brighter that 22.2 mag (the 95% completeness 
limit of the photome tric survey) were co-a dded. We note 
that modeling of the iBlanton et al. 1 (I2003T) SDSS r band 
luminosity function indicates that even in our most dis¬ 
tant clusters the contribution to the total light of galaxies 
fainter than our upper limit of 22.2 would be less than 2.5% 
(i.e a log(luminosity)=0.01). This contribution is therefore 
deemed insignificant to our results, and is ignored. 

A colour cut excluding all galaxies 0.2 mag redder in 
(g — r) than the BCG was applied in order to minimize 
the impact of contamination by high redshift interlopers. 
In order to estimate the background contamination, the 
same process was applied to seven concentric annuli between 
3R200 and 4i?2oo • These annuli were selected to have an area 
equal to the area of the central region within R. 200 - In each 
of the g, r and i bands the median value of the seven annuli 
was then taken as the background level and this value was 
subtracted from the value calculated within R200 yielding 
an estimate of the total luminosity of faint galaxies within 
77200 - These luminosities were a lso then k-corrected to zer o 
redshift using the functions of IChilingarian et al. I (I2OI0I) . 
The total luminosity in faint galaxies was then added to the 
total luminosity in spectroscopically confirmed galaxies to 
yield a total luminosity for each group. 

However, in some cases the background luminosity in 
faint galaxies exceeded the estimated total luminosity of 
all the galaxies in the cluster, resulting in negative values 
for the luminosity of the cluster, clearly indicating that 
the background level was overestimated. Such clusters were 
excluded from our analysis. In other cases, despite high 
background levels, the total luminosities of the clusters 
remained positive raising the possibility that the clusters 
had, indeed, very low luminosity compared to the back¬ 
ground against which they are projected. In such cases 
we excluded only clusters in which the background level 
exceeded five times the luminosity of faint galaxies in the 
cluster, deeming the uncertainties in the background level 
to be too large for a reliable total luminosity estimate. 

After all the selection cuts above were made we 
retained a total of 614 clusters. 

In order to compare literature values in the B and 
K bands we used the colour depen dent functi o ns o f 
IBlanton fe Roweis I (120071) (B band) and lYaz et al. I (l20ld) 
(K band) to estimate the luminosities of each galaxy in each 
cluster in these bands. The total luminosities in the bands 
were then also estimated as above. 
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3 RESULTS 


In this section we report our results and compare them 
to five recent works from the literature that measure 
masses and luminosity ratios in significant samples of 
groups/clusters in the B, r, i and I\ bands (recall that, 
for our results, r and i are taken directly from the SDSS 
database while B and K are generated from the SD SS data 
using the c olour dependent t ransforms of lBlanton fe Roweis I 
(l2007f) and lYaz et al. (H) respectively). A full table of re¬ 
sults is available from the authors and will be made available 
on-line. 

Fig [3] shows mass against luminosity for the clusters 
in our final sample in all four bands (436, 614, 611 and 
427 clusters in B, r, i and K bands, respectively). Lines of 
constant log of the mass-to-light ratio (\M/L\=\og(M/M q)- 
log (L/Lq)) are shown for [M/L]=3.0, 2.5, 2.0 and 1.5 (de¬ 
creasing to the lower right). 

In the figure our data are shown as black points, while 
in each panel, data from the literature are also presented 
as either red points (when actual data are available) or red 
lines (when the literature provide fits to data). Note that 
in all cases the literature data have been converted (when 
necessary) to a Ho=70 km s -1 Mpc~ 4 , S2a=0.7 cosmology, 
k-corrected to z = 0 and adjusted the solar luminosities 
given in the introduction. 

Fig. [3] generally shows good agreement between our 
data and the literature studies shown in red, particularly 
at masses higher than >10 14 Mq. However, there appear to 
be two areas of poo rer agreement - the fi rst in the compar¬ 
ison to the study of iGirardi et al. 1 (|2002h . and the other at 
the low mass end of our mass range (i.e. <10 14 Mq) in the 
B , r and i bands. 

Addressin g first the poor agreement with the study of 
IGirardi et al. 1 120021 in the B band , there is clearly an of fset 
between our data and the results of lGirardi et al. 1 1120021 ) at 
higher masses (>10 14 Mq), with the average of log(L) within 
±0. 1 dex of log(M ) =14 .25 being 11.67 and 11.97 for our 
and IGirardi et al. I l2002h data respectively - i.e. a 0.3 dex 
(or factor of 2) difference in luminosity at a given mass, and 
therefore in the mass-to-light ratio. However, we note that 
our results are consistent with the study of E ke et al (2004b ) 
in the B band, as well as wit h the studies of Popesso et al. I 
d2007tl , I Sheldon et al. I (120091 4 , iRamella et al. I J2004I ') in the 
r, i and K bands respectively. We then conclude that the 
Girardi et al. (2002) study suffers from an, as yet unidenti¬ 
fied, systematic difference in their luminosity estimates with 
respect to other recent studies. We assume the offset to be 
in the luminosity direction because, as noted in Section 12.41 
we have already shown that the mass estimation method¬ 
ology of lGirardi et al. I d2002l ') is completely consistent with 
the simple dynamical mass estimates made in this paper. 

Turning our attention to the disparities at low masses, 
our results clearly show a pronounced increase in the slope 
of the mass-luminosity relation as we go to low masses. How¬ 
ever, we note that many of the literature studies presented 
in Fig. [3] appear to show similar, albeit smaller, increases 


in the slope of the mass-luminosity relation a t the low mass 
end. The effect is c l early se en in the studies of lGirardi et al. I 
(20021 ,lEke et al~l ( 2004b|) and can also be seen in the data 
of ISheldon et al~T ( 20091 1 in which the gradient decreases 
from ~1.6 below 10 14 Mq to ~1 . 3 abo ve this mass. The 
power -law fits of of lPopesso et al.l (120071) or IRamella et al!~l 
(|2004h do not, per force, e xhibi t su ch gradient changes. How¬ 
ever, in the case of lRamella et al.~l ll2004lf . the study contains 
ver y few points at the l ow masses of interest here, while 
the lPopesso et al~l J2007h study presents data selected from 
two catalogues, one opti cal, the other a n X-ra y catalogue. 
Examination of Fig. 9 of lPopesso et al.l (|2007l ) shows that 
the optically selected sample does exhibit an increase in the 
slope of the mass-luminosity relation at the lowest masses, 
while the X-ray selected systems do not. This suggests that 
this trend may be, at least in part, the result of optical selec- 
tio n effects (a p oint that was also raised in the recent paper 
of iMulrov et al. I d2014|j 'l. In fact, our study, which requires 
at least 6 bight galaxies within R 200 of galaxies in the red- 
shift range z = 0.05 to 2 = 0.16, is generally biased towards 
higher redshifts than the other studies which are dominated 
by much closer systems. This raises the possibility that this 
optical selection effect may be playing a stronger role shap¬ 
ing the locus of our results than in the other studies pre¬ 
sented here. To examine this possibility in more detail we 
again turn to the Millennium Simulation. 

In Fig. [4] we show the mass-luminosity relation for 
groups and clusters in the Millennium Simulation. The me¬ 
dians of the dynamical mass data in bins of luminosity are 
shown as solid red lines. The 15th and 85th percentiles are 
shown as dashed lines. In the top panel we plot the masses 
and luminosities taken directly from the simulation. In the 
bottom two panels we plot dynamical masses derived from 
the simulation (see Section [22]). In the top and middle plots 
all groups and clusters with 6 or more member galaxies (ir¬ 
respective of their luminosities or redshifts) are plotted. In 
the bottom panel we apply the SDSS selection criteria - i.e. 
6 or more galaxies brighter than 17.7 mag in the r band, 
and redshifts in the range 0.05< z <0.16. The bottom panel 
therefore represents the Millennium Simulation prediction 
for the locus of our observational data. 

Comparing first the top two panels, the displacement 
and non-linearity in the relation between Millennium Simu¬ 
lation mass and dynamical mass identified in Fig. [2] is also 
evident when comparing these two plots. A marked increase 
in scatter is also clearly seen in the dynamical mass data. 

Comparison of the bottom two panels shows the effect of 
requiring 6 or more bright members (r < 17.7 mag) within 
R 200 and restricting the redshift range to 0.05< z <0.16. 
Clearly, at low masses, when combined with the require¬ 
ment of a minimum number of bright members for dynami¬ 
cal analysis, this redshift restriction causes only the brightest 
groups and clusters at any given mass to be selected, thus 
introducing another strong non-linearity into the apparent 
mass-luminosity relation. What is more, comparison of the 
bottom plot with the r band data in Fig. [3] suggests that 
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log(Ltot) (L©) logCl^J (Uq) 


Figure 3. Mass is plotted against total group/cluster luminosity. Black dots are our data, while red points and lines are taken from 
literature studies (as detailed in the bottom right of each panel). The four black lines in each panel are lines of constant log of the 
mass-to-light ratio([M/L]), with values of 3.0, 2.5, 2.0 and 1.5, decreasing from upper left to the lower right. Agreeme nt between our 
data and literature data is generally good at high masses, but at lower masses and in the comparison to lGirardi et al.~l l2002i 'l in the B 
band discrepancies are evident. 


Table 1. The offsets from the median of the 15th and 85th per¬ 
centiles in three masses bins in both our data and the Millennium 
Simulation (MS). 



Data 

MS 

log(M) (M 0 ) 

15% 

85% 

15% 

85% 

13.3 

-0.52 

0.39 

-0.50 

0.45 

14.0 

-0.32 

0.30 

-0.39 

0.33 

14.7 

-0.27 

0.27 

-0.28 

0.25 


the Millennium Simulation is reproducing the observational 
data extremely well. 

In order to demonstrate this, in Fig. [5] we plot the pre¬ 
dictions of the simulation for the mass-to-light ratio against 
mass scaling relation and compare to the medians of the ob¬ 
servational data and their uncertainties. We choose to show 


mass-to-light ratio here as the relatively small dynamical 
range of this parameter serves to emphasize similarities and 
differences between the various data sets, while the binning 
of the data has significantly reduced the correlated errors 
which would otherwise make interpretation of such a plot 
difficult. 

In Fig. [5] the data taken directly from the Millennium 
Simulation is shown as a thick solid line. For ease of com¬ 
parison, this line has been shifted by the ~0.1 dex offset 
between Millennium Simulation and dynamical masses first 
identified in Fig. [2] and also evident in Fig. [4] The data for 
dynamical masses with no luminosity or redshift selection 
criteria applied is shown as a thin solid line, while the dy¬ 
namical mass data with luminosity and redshift selection 
criteria applied is shown as a dashed line. 


© 0000 RAS, MNRAS 000, 000-000 
















8 Proctor, Mendes de Oliveira, Azanha, Dupke, Overzier 





Figure 4. The r band mass-luminosity relation of groups/clusters 
in the Millennium Simulation. In the top panel Millennium Sim¬ 
ulation mass (rather than dynamical mass) is plotted, while the 
bottom two panels show dynamical mass. In the top two panels 
all systems with redshift less than 0.16 are shown, while in the 
bottom panel our observational selection criterion of a minimum 
of 6 bright galaxies within R 200 of systems at 0.05< z <0.16 is 
applied. Solid and dashed red lines are the median with 15th and 
85th percentiles, respectively. 



Figure 5. Mass-to-light ratio in the r band are plotted against 
against mass. The medians of observational data and their (cor¬ 
related) uncertainties are shown as black points with error bars. 
Millennium Simulation data are shown as red lines. Data taken 
directly from the Millennium Simulation are shown a thick solid 
line, while the data for dynamical masses derived from the simu¬ 
lation with no selection criteria applied (see text) are shown a a 
thin red line. The data from the simulation using dynamical mass 
and with selection criteria applied is shown as a dashed line. The 
simulation clearly reproduces the observational data extremely 
well. 

What is more, in Table[l]we show the width of the 15th and 
85th percentiles in dynamical mass at selected points along 
the mass-luminosity relations of the simulated and real data, 
and agreement is again excellent. Clearly, the Millennium 
Simulation is doing a remarkably good job of reproducing 
the actual observational data. 

Our analysis of the Millennium Simulation has therefore 
clearly identified the importance of selection effects in shap¬ 
ing the locus of the mass-luminosity relation found from our 
data, and suggests that such effects may be present to, some 
lesser extent, in the other studies considered here. Indeed, 
the excellent agreement found between our data and the 
Millennium Simulation suggest that the simulation might 
be used to identify, and even correct for, such effects in fu¬ 
ture studies. 


The changing slope in the Millennium Simulation mass 
against dynamical masses relation noted in Fig. [2] at ~10 14 
M 0 is again evident in this plot, with the two solid lines ex¬ 
hibiting the same shallow slope above this mass, but steep¬ 
ening and diverging below. Examination of Fig.[5]also shows 
that the Millennium Simulation prediction for the locus of 
mass-to-light ratio against mass relation for data using dy¬ 
namical masses and with selection criteria applied (dashed 
line) is an extremely good match to the observational data. 


4 CONCLUSIONS 

As the first step in a program to perform a detailed 
study of “fossil groups”, we have estimated masses and to¬ 
tal luminosities in ~600 groups and clusters identified in 
the maxBCG catalogue using SDSS photometry and spec¬ 
troscopy. We have also performed an analysis of Millennium 
Simulation data subjected to the same processes and selec¬ 
tion criteria as the observational data. We have compared 
our results with five other studies in the literature in four 
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different wavebands ( B , r, i and K) and the Millennium 
Simulation data. 

At high masses (>10 14 Mq) we find extremely good 
agreement with four of the five literature and the Millen¬ 
nium Simulation results, but find an exception in a single 
literature study llGirardi et al. 1 12002| ). which appears dis¬ 
cordant. At lower masses the agreement of our data with 
literature results is poorer. However, we have used the Mil¬ 
lennium Simulation to demonstrate that a large part of the 
discrepancies at these masses can be explained by selection 
effects. 

We have also used the Millennium Simulation to pro¬ 
vide a hint to the cause of remaining discrepancies at 
lower mass, showing that within the simulation the as¬ 
sumption that mass is proportional to cr 3 begins to fail as 
one approaches lower masses. With hints of a similar phe- 
nomena already hav ing been observed in X-ray data (e.g. 
iHelsdon & Ponman Il2000l ; lOsmond & Ponman lEpol l. this 
is an issue that certainly warrants further investigation. It 
should be noted that our findings with regard the Millen¬ 
nium Simulation are of particular importance in studies that 
attempt to measure the gradient of the mass-to-light ratio 
against mass scaling relation which, as is evident from our 
results, are extremely susceptible to these effects. 
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